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GLOBAL STABILITY OF STEADY TRANSONIC EULER SHOCKS IN 
QUASI-ONE-DIMENSIONAL NOZZLES 

JEFFREY RAUCH, CHUNJING XIE, AND ZHOUPING XIN 

Abstract. We prove global in time dynamical stability of steady transonic shock so- 
lutions in divergent quasi-one-dimensional nozzles. We assume neither the smallness of 
the relative slope of the nozzle nor the weakness of the shock. Key ingredients of the 
proof are an exponentially decaying energy estimate for a linearized problem together 
Xv ' with methods from IT21 . 

< 

Compressible isentropic Euler flows in quasi-one- dimensional nozzles are governed by 

. a'(x) 
Pt + [pu) x = r-/«, 

a'x ■ (1) 

\Q ■ (pu)t + (pu 2 + v{p))x = — P« 2 

in ' 

oo 
in 



o 



£ 



1. Introduction and Main Results 



a(x) 
where p, u and p denote the density, velocity and pressure, and, a(x) is the cross-sectional 
j--. area of the nozzle. The typical examples are ID, 2D rotationally symmetric, and 3D 



spherical symmetric compressible Euler equations, where a(x) are 1, x and x 2 , respectively. 
We assume that p satisfies 



Vp>0, P (p)>0, p'(p)>0, and, p"(p)>0. 

Examples are p(p) := Ap 1 with 7 > 1, in particular the polytropic (7 > 1) and isothermal 
(7 = 1) gases. The sound speed is c(p) := \/p'{p)- A state p, u at which the flow speed 
is larger than the sound speed, \u\ > \Zp'(p), is called supersonic. The opposite case, 
M < \/p'{p)) is subsonic. 

For steady states, (pj) becomes the system of ordinary differential equations 

a'(x) 
(pu) x = r^P«, 

<M • < 2 > 

{pu 2 + p(p)) x = r^pu 2 ■ 

a(x) 

studied intensively in [HE]. Some properties needed later are the following. 
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The first equation in (J5J), (apu) x = 0, holds if and only if 

3M 6 R, apu = M . (3) 

After a possible reflection we may assume without loss of generality that M > 0. The 
flow is from left to right. Using ([3]) to eliminate u from the second equation in ([2]) yields 

a scalar equation for p, 

fM 2 \ 

+ ap(p) x = 0. (4) 

V a P J X 

Definition 1. A steady transonic shock solution is a piecewise smooth solution of Op 

with two smooth solutions separated by a shock connecting a supersonic state on the left 

to a subsonic state on the right. 

The steady transonic shock solutions of ([2]) on [l,L] are the piecewise smooth solutions 

_ f (p_, M/(ap-)), HI <x <x , 
\ (p + ,M/(ap + )), iix <x<L, 

satisfying the Rankine-Hugoniot condition 

P(P-) + / - p j (»0~) = \P(P+) + 7— -p J ( X + ), (6) 

with p_ is supersonic (p'(p-) < , ^ )2 ) and p + is subsonic (p'(p + ) < ,^1 )2 ). This implies 
the entropy condition 

p-(x ) < p+(x ). (7) 

Steady transonic shock solutions were constructed in [Tj, [U] satisfying the boundary 

conditions 

M 
(p,u)(l) = {pi,u{) = {pi, -tjt— ), and, p(L) = p r , (8) 

o(()pi 

with (pi,Ui) and (p r ,u r )) supersonic and subsonic respectively There have been many 
studies of the stability of steady transonic shocks for flTJ (see [ID])- In |TUJ a wave front 
tracking variant of Glimm's scheme was used to prove that when \a'{x)/a{x)\ is small, 
a weak transonic shock is dynamically stable if a'(xo) > and dynamically unstable if 
a'(xo) < 0. Studies on the solutions of general hyperbolic conservation laws with source 
terms include [HI [HI El El H] and references therein. For piecewise smooth initial data, in 
the case that a(x) = x or a(x) = x 2 , Xin and Yin [18] proved dynamical stability of weak 
shocks in nozzles with large / (thus \a' (x) / a(x\ small) in the class of piecewise smooth 
transonic shock solutions. The main result of this paper proves global stability in time 
without these restrictions. 



Let (p, u) be a steady transonic shock solution of the form (jSJ) satisfying the boundary 
conditions (jHJ). Suppose that the solution stays away from vacuum, 

inf p(x) > 0. (9) 

xe[i,L] 

Solving an initial value problem for fl2]) we can extend (p_, uJ) to be a smooth supersonic 
solution to (J2J) on [l,Xo + S] for some S > 0, that coincides with (p_,w_) on [Z,£o]- The 
notation (p_,u_) is used for the extended solution as well. Similarly, denote by (p + ,u + ) 
a subsonic solution of (J2]) on [x — 5, L] for some 5 > 0, that coincides with (p + ,w + ) on 

[x ,L]. 

Consider the mixed initial boundary value problem defined by ([I]) with initial data, 

(p,u)(0,x) = (po,Uo)(x), (10) 

and boundary conditions 

M 
(p,u)(t,l) = KPh-^Wz)i P^ L ) = Pr> (11) 

with pi and p r from (jHJ). 

Assume that the initial data satisfy 

{(p -,Mn-)(x), if I < X < Xn, 
(12) 
(po + ,u +)(x), if x < x < L, 

These data are small perturbations of (p, u) in the sense that 

\xq-x q \ + \\(p 0+ ,u 0+ ) - (p+,u+)\\ Hh+a( p 0tL]) 

+ \\(po-,uo-) - (p_,M_)|| Hfc+2(Mo]) < e, 

for some small e > 0, and some integer fc > 15, where xq = min{xo,xo} and xq = 
max{a;o,5o}- Moreover, (pq,Uq) is assumed to satisfy the Rankine-Hogoniot condition. 
That is, at the shock location x = Xq, 

{(p(Po+) + Po+u 2 0+ - (p(po-) + Po-«o-)) (Po+ - Po-) = (Po+%+ - Po-%-) 2 - (14) 
Our dynamical stability theorem is the following. 

Theorem 2. Let (p,u) be a steady transonic shock solution to system |lj) satisfying |3]) ; 
(Eli, $E), and (G|). TTie &ey hypothesis is that the nozzle is widening at the shock location, 

a'(x ) > 0. (15) 



(13) 
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There exists an e > so that for any < e < e , if the initial data (p , u ) satisfy ip~3\) , 
P7P an d the compatibility conditions at x = I , x = xq and x = L or order (k + 2), then the 
initial boundary value problem (QP, / fiOj) and < f77]j has a unique piecewise smooth solution 
(p,u)(x,t) for (x,t) e [l,L] x [0, oo) containing a single transonic shock x = s(t) with 
s(0) = xq and I < s(t) < L satisfying the pair of Rankin-Hougoniot conditions for t > 0, 

(p{p) + pu 2 ){t, s{t+)) - {p{p) + pu 2 ){t, s{t-)) = (pu(t, s{t+)) - pu{t, s{t-)))s{t), 

pu(t,s(t+))-pu(t,s(t-)) = (p{t,s(t+))-p(t,s(t-)))s(t), 
and the Lax geometric shock conditions, 

(«-v9G?)(*, *(H) > Kt) > (u-^{pj)(t,s(t+)), (u+^p 7 (p))(t,s(t+)) > s(t). 

Denote 

] (p-,u-), if I < x < s(t), 
{p,u) = < 

{ (p+, u +), if s(t) <x< L. 

There is a To > so that for t > Tq, 

(p-,U-)(t,x) = (p-,U-)(x), for I < x < s(t) . 

The solution approaches the steady transonic flow at an exponential rate, that is, there 
exist positive constants C > and A > independent of the solution so that for t > 

fc-6 

\\(p + ,u+)(;t) - (p + ,u + )(-)\\ wk -7,oo {s{t)>L) + J2 \dT(s(t) - x )\ < Cee- X \ 

m=0 

where {p±,u±) is the unperturbed solution. 

Remark 1. The results in Theorem^ are also true if we impose small perturbations in 
the boundary conditions |2j). 

Remark 2. The regularity assumption in |73|) is not optimal. Adapting the paradifferen- 
tial methods from [H], one can decrease the regularity required in |73j). 

Remark 3. Compatibility conditions for the initial boundary value problems for hyperbolic 
equations are discussed in detail in [T6| 



Remark 4. We require neither the smallness of \a'/a\, largeness of I, nor the weakness 
of the shock strength as in [18J . Similar results for steady transonic shock solutions for 
the Euler-Poisson equations were proved in [12] . The strategy of the proof of Theorem 1 
is inspired by [T2] . 



Remark 5. If the condition < f73]) is violated, i.e. o'(xq) < 0, then the transonic shock 
solution is unstable, see [TU1 fl8] . 

The paper is organized as follows. Section 2 transforms the problem to free boundary 
problem for a second order scalar hyperbolic partial differential equation. The weighted 
energy equality for the associated linearized problem is derived in §3. This yields Theorem 
1 with the aid of ideas from 1121. 



2. Transformation of the Problem 

Let (p,u) be a steady transonic shock solution of the form (T5]) satisfying (EJ), (IE]), OH]), 
and (J9]). If the initial data (po, ^o) satisfies f|T3]) and the compatibility condition of order 
k + 2. It follows from the argument in [7] that there exists a piecewise smooth solution 
of the equations on [0, T] for some T > 0, which can be written as 

k u)= r(p-,«_), if i <*<*(*), (16) 

I (p+, w + ), if s(£) < x < L. 

Since the flow is supersonic on the left it follows that there is a To > so that, when 
t > T , (p_, uJ) depends only on the boundary conditions at x — I. Moreover, when e is 
small, by the standard lifespan argument, we have T < T (cf. [7]). Therefore, 

(p-,u_) = {p-,uJ) for t > T . (17) 

It suffices to study the solution for T > so without loss of generality we suppose that 
To = 0. We want to extend the local-in-time solution to all t > 0. In view of ffTTl) . it 
suffices to obtain uniform estimates in the region x > s(t), t > 0. We formulate an initial 
boundary value problem in this region. First, the Rankine-Hugoniot conditions for (TT6l) 
read 

[pu] = [p]s'(t), [pu 2 +p] = [pu]s'(t), (18) 

using standard notation for the jump [/] :=, f(s(t+),t) — f(s(t—),t), so 

[p + pu 2 ] [p] = [pu] 2 . 



Denote M := apu. Then 

/ w , ^ Ml{t,s{t)) . ..... Mi 

p(p+)(t,s(t)) + -—±^-ll—- p (p + )(s(t)) - -—±—( s (t)) 



a 2 (s(t))p + (s(t) : t) 



a 2 (s(t))p^ 



M 2 M 2 

+ p(/*f )(*(*)) + ^W*)) " P(P-Mt)) - — (*(*)) ) (p+& *(*)) - p-(s(t))) 



a?p + 



a?p- 



M + (t,s(t))-M+(s(t)) 



a(s(t)) 
The implicit function theorem and the momentum equation in (J2J) imply, 

(M+ - M+)(t, s(t)) = M((p+ - P+)(t, *(*)), s(t) - x ) (19) 

where ^i regarded as a function of two variables satisfies ^/\(0, 0) = and 



dsrfi 



d{p+ - p H 



a(p'(p+) ~ u 



r,2 



( x o) 



9M 



(0,0) 



o'w_(p + -p_) 



(z )- 



(0,0) 



2ii4- 9(s — xo) 

Note that M + = M_ = M. And so substituting ([19]) into the first equation in ([TBI) yields, 

s'{t) = s6(p + -p + ,8(t)-x ) (20) 

where s^ 2 satisfies ^(0, 0) = and 



dss/o 



d(p+ ~ P+) 
Define 



r,2 



(0,0) 



P'(p+) - u 
2u+(p+ - p_) 



ft^o 



■^0 



d(s(t) - x ) 



a'u. 



(0,0) 



2a 



^o 



?/;_(£, x) : = -Mt+ J a(y)p-(y)dy 



and 



4>+(t,x) := ip-(t,x )+ I a(y)p + (y)dy, ip+{t,x) = ip-(t, s(t)) + / a(y)p + (y,t)dy. 

J 'if) 



Set 



Then 



* := ^ + (x,£)-^+(£,x). 



* t = M - M+, V x = ap + - ap + . 
The second equation in ((Tj) is equivalent to 

-(V>+ + V) tt + 9 S — — = — + a(x)d x p I — I = 0. 



<9^+ + 9** 



ax 



Write this equation as follows. Setting £ = (£o>£i) = (t,x), 

where a^-, 6, and c are smooth with respect to each variable, and satisfy (using the Einstein 
summation convention), 

^ m = ciijix, 0, 0)8^ + bi{x, 0, 0)<9;# + cg{x, 0, 0)* 

= d tt ^> + 2u+d tx V + 04 - p'(p + ))d xx V + 2d x u+d t V + Bd x y 

with 

B ■= -p'(p + )-d x p'(p + ) + d x u 2 + . 
a 

In terms of \1/, the equations ( 1T9|) and ( 120]) can be written as 

%(t, s(t)) = -M( *f' s( * )} , S (t) - x ), (22) 

and 

^ = ^(^7?P,*(*)-Xo), (23) 

respectively. A direct computation yields 

*(t, s(t)) = </>+(£, s (t)) - $+(*, s(t)) = ^_(t, s(t)) - Mt, s(t)) 

= ip-(t,x ) +ip-(t,s(t)) -ip_(t,x ) -tp+(t,x ) -ip + (t,s(t)) +$ + (t,x ) 
= d x ip_(t,x ) ■ (s{t) -x ) -d x ip + (t,xo) ■ (s(t) -x ) +R U 
where R\ is quadratic in s(t) — xq. This implies 

s{t)-x = tf z {V(s(t),t)), (24) 

where srf-& satisfies ^(0) = and 



It follows from (l22]h (123]) and (124) that 



o(P- - P+) 



(x ) 



d t m = £/ 4 (^ x , *), at x = s(£), (25) 

where 

^ 4(0 ' 0) = ' 0*1 = 2, + (Xo) ' -W = -^ {Xo) - 

On the right boundary, x = L, ^ satisfies 

d x ^> = 0, at x = L. (26) 



It suffices to derive uniform estimates for \l/ and s that satisfy fl2Tl) . fl24|) . ( 125]) and ( 126|) . 
To this end, first transform this free boundary value problem into a fixed boundary 
value problem. Set 



t :=£, 



, T s X — s(t) 
X := [L — So) - ? 77T + x 



L - s(t) 



and 



Then 



qi{x,a) 



L — x 



L — xq — a(t) 
d d . ~ d 

m = dt ~ ° {t)qi Yx 



g 2 (<x 



a(i) := s(t) -x , 

L — Xq 

L — Xq — a(i) 
d d 

■5- — <?2-^r, 
OX ox 



(27) 



d 



w = w + ^w* - 2 ^®4- - * O + 2 l^S[ f) ; * 



J 2 



dxdt 
So (12TT) becomes 



?2 



/ <9 2 



</(t) g d 2 

V<9x<9£ + L - xq - a{i) dx qx ° W dx 2 



8F_ 2 d?_ 

~ ?2 Xz.2 ■ 



%* + 92% ■ — t + a{x)q 2 d i p 

ap + + q 2 d^ 



dx 2 dx 

ap + + q 2 ^x 



a[x) 



~ 2a'{t)q l d^ + ( qi a'(t)) 2 d^ - 2 t ^'^ .-. q&V = a"(i) qi d^. 

L — Xq — a[t) 

Equation (I24|) takes the form 

<T = */ 3 (V(t,X = X )), 

and the equation for the shock front, ( 123|) . becomes 

'q 2 (ff)V 



da 
dt 



a(t) 



Using (liMl) to represent the quadratic terms for a in terms of \1/, we have, at x — Xq, 

da a'u. 



where ^2 satisfies 



^ 2 (*s,*) + 



dt 2a 

c 2 (p + ) 



(X y = «2(*5,*), 



(2f 



M 



-(^o)*d 



<C(^ 2 + ^ 2 ). 



2 (p+ - p-)u+a 

It follows from (12^1) and ( 1281) that one can represent a and a' in terms of \I/ and its 
derivatives at x = xo- Thus, after manipulating (|25|) with ( ]2"%|) and (|28|) . one finds 

% = «i(*2,*), at £ = x . 



By the implicit function theorem this is equivalent to 

*« = ^3(*f,*), at x = x , 

where ^3 satisfies 



%(%,*) 



2u_ 



M*t 



a m_Mj 



(x )# 



< C(^ + ^ 2 ). 



c 2 (p+) - u\ v " u; * (c 2 (p+) - u%)a ' 

In the following, drop the in x and t for ease of reading. 

In summary, the problem has been transformed to the following compact form 

' &{x,y,a)-y = a"{t) qi d x ty, {t,x)e[0,oo)x[x ,L}, 

d x ^ = rf 1 (^ t ,^)^ i + e 1 (^ i ,^)^, atx = x , 

d x ^ = 0, at x — L, 

o(t) = s/ 3 (*(t,X )), v{0) = a , 



where, 



i,j=0 



j=0 



with 



di(*t,* 



1 && 



Furthermore, one has 



(e%,e*)de, e x (**,*) 



Sf(x, 0, o)$ = J&f $ 



^ 



(e%,ev)de. 



and 



a o(a;, V\l/, a, cr') = 1, Ooi(^) 0, 0, 0, 0) = a w (x, 0, 0, 0, 0) = w H 

a n (x,0,0,0,0) = -(p'(p + )-^), 

b {x, 0, 0, 0, 0) = 2d x u+, h(x, 0, 0, 0, 0) = B, 

1u+ a'(xo)u+U- 



di(0,0) 



l^oj 



ei(0,0) 



c 2 (p+)-w+ ' '"" (c 2 (p + )-u 2 + )a 

3. Proof of the Main Theorem 



[xq). 



(29) 



The key first step in the proof of Theorem [2] is the exponential decay of solutions of 
a linearized problem. It yields in non trivial manner a priori estimates for the nonlinear 
problem which lead to the global stability of the steady shock. 
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3.1. Linear Estimate. The first step is the following energy identity for the linearization 
at the steady shock. 

Lemma 3. Assume that a satisfies < f73j) . Let ^ be a smooth solution of the linearized 
problem 



' Jzf ^ = 0, x < x < L, t > 0, 



d x v 



2ua 



c 2 (p+ 
cL* = 0, at x = L 



. a'(xo)u + U- 



x )ty, at x — xq, 



^(0,x) = hi(x), *t(0,x) = /^(z), x < x < L. 



(30) 



T/ien t/ie following dissipation identity holds, 



E(*,t) = E{*,0)-D{%t), 



(31) 



where E and D are defined as, 



E(*,t) 



a'uiu-^ 2 



(t,x )+ / u + {(d t ^) 2 + (p'(p + )-u 2 + )(d x ^) 2 }(t,x)dx. 



■I'D 



DOM := 2/ t£(xo)(W(t,z o ) + <(20(W(t,£) dr. 



Proof: Multiplying the first equation in ( I30p by ■u + (;r)c? t \I/ and integrating by parts yields 



t rL 



&o$ ■ M+9 t ^ dxdr 



./no 



w + (^) 2 H-«+(c 2 (p + -^)( 9 ^)^ 



(£, a:) dx 






/ / [Bu + + 9 x (M + (c 2 (p + ) - ul)))d t Vd x y dxdr 

J0 Jx 



' - s x=L 

T-.2 fO 1 Tr^2 i „-; fr-,2 "2 { - 



(32) 



+ / (u 2 + (d t yy + u + (u 2 + -c 2 (p + ))d x m t y 



dr 



X=XQ 



L „-; 



u + (d t vY + u + (c 2 (p + -ui)(d x *y 



(0, x) dx . 



.1-0 
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Note that 



Bu+ + d x (u + (x)(c 2 (p + ) - u\)) 

= Bu+ + u + (x)d x (c 2 - uV) + d x u + (x)(c 2 - uV) 



U A 



- P (p + )-p(p + )^ + H ^j 



+ u + 9 K U'(p H 



(33) 



M 2 



(0*Vm 



+ d x u + (p(p) -u\) 



a ax \ap^ ' 



The momentum equation can be written as 



Or 



M 2 
a{x)p. 



+ a(x)p'(p+)p' + = 



so, 



P'(p-\ 



M 2 



(ap + ) 2 



P+ 



M 2 
a 3 p + 



a'. 



Therefore, 



The boundary terms in (1321) are 



a' ,, „ d ( M \ , 
u+-p'(p + ) + — — * ' ' - > ^ 
a ax \ a Pn 



(p'(p + )-4) = 0. 



(34) 



uKx)^*) 1 



x=L 



X=XQ 



+ u+(ui - c 2 )d x mv 



x=L 



X=XQ 



= U 2 + (d t ^Y(T,L) +ui(x )(d t ^Y(T,Xo) + 



a' u\_U- 



a 2 



■d t ^ 2 (r,x ) 



(35) 



= u 2 + (d t ^) 2 (r, L) + u 2 + (x )(d t ^) 2 (T, x ) + d t 
The lemma follows from (j32|), (|33|), (134J), and (135|). 



2a 



(r,x ) 



a 



The estimate in Lemma [3] and the method of [12] imply that solutions of the linear 
problem ( )30l) decay exponentially to zero. This in turn implies the exponential decay of 
solutions of the nonlinear problem. We describe the form that these arguments take in 
the present context. 

Lemma 4. If a' satisfies < f73|) then there exist constants Xq > and C > so that 
solutions \& of the problem §W\) satisfy, 



E(V,t) < Ce- Xot E(ty,0) 
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Proof: Step 1. Rauch- Taylor estimates. Thanks to the boundary condition at x — Xq, 
the estimate (13TT) implies 

E(*,t) H-Ci [ (^ 2 t +^l)(s,x )ds < E(^,0) + C 2 [ ^ 2 (s,x )ds. 
Jo Jo 

The argument in [T7] and the details in [TJl Appendix], one can show that there exist 

T > and 5 e (0, T/4) such that 

f (V 2 t + *l)(t,x )dt >5E(^,T)-C 3 [ V 2 (t,x )dt. (36) 

Jo Jo 

Combining ( 13T1) and ( 1361) yields, 

(H-C 4 )£(*,T) < E(V,0) + C 5 [ V 2 (t,x )dt, (37) 

for some positive constants C 4 and C 5 , independent of t. 

Step 2. Spectrum of the evolution operator. Define a new norm || • ||x for the function 
h = (h x ,h 2 ) CH 1 xL 2 ([x ,L]), 

\\h\\x = ^ L ^M\hi\ 2 M+ / u + {\h 2 \ 2 + (p'(p + )-u 2 + )\h[\ 2 }(x)dx. 

a Jx 

The associated complex Hilbert space will be denoted by (X, || • || x ). Define the solution 
operator St '■ X *-$■ X as 

St(h) = (9(t,-),%{t,-)) 

where \1/ is the solution of the problem (1301) with the initial data h = (hi, h 2 ). Applying 
the Lemma on page 81 in |15| . there are at most a finite set of generalized eigenvalues for 
the operator St in the annulus {jtc~ < M < 1} C C, each with finite multiplicity. 

Step 3. Refined estimate for the spectrum of St- We show that the spectrum does not 
touch the unit circle. Otherwise there would exist wel and VeX such that 

(S T ~ e iuJ I)V = 0. 

Note that the identity (I3TI) still holds in the complex setting if we replace the square terms 
in E and D by the square of modulus. Thus 

E(V, 0) - E(V, nT) = nD(V, T). 

Since E(^>,nT) and E(^>, 0) are both positive and finite, it follows that 

£>(tf,T) = 0. (38) 
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Therefore E(V, t) = E(V, 0), for all t. Let V := ker(S T - e iuj I). Note that the coefficients 
in the problem (1301) do not depend on t, so 

{S T ~ e w I)S t = S t (S T - e iu I) 

In particular, S(t)Y C V for any t, so V is invariant with respect to St. Therefore S t \\ 
is a semigroup on a finite dimensional subspace. This yields that 5"(t)|v = e tA for some 
A G Horn V. The definition of V implies that e T A = e tuJt I so the spectrum of A is purely 
imaginary. Choose an eigenvector w of A, Aw = Xw, X G i(3 G iR Then S t w = e xt w. 
From Stw = e tl3t w and (|38l) it follows that w(xo) = 0. The boundary condition at xq 
then yields w'(xq) = 0. The uniqueness of solutions of the linear homogeneous second 
order ordinary differential equation satisfied by eigenfunctions implies that w — 0. This 
contradicts the assumption that w is an eigenvector hence not equal to zero. 

The contradiction shows that there exists < f3 < 1 so that the spectrum, <j(St), 
of St satisfies ct(St) C {|z| < /3 }. The formula for the spectral radius implies that 
\\S(nT) || = HS^T)™!! decays exponentially as n — )■ oo. This is equivalent to the assertion 
Lemma HI □ 

Corollary 5. Assume that a' satisfies l[T5\) and < k G N. Define 

k 
E k {*,t) := ^2E(8F%t). 

m=0 

Then, with Xq from Lemma\l^ solutions \l/ of the linearized problem (|30|) satisfy 

E k (V,t) < Ce- Xot E k {V,0), and, 

/ e Ao4 / 4 ^(|9^| 2 (t,a;o) + |9^| 2 (t,L)) dt < C^(f,0). 

Proof: Since the coefficients of the equation and the boundary conditions are independent 
of t, applying this estimate to the solution <9™\l/ yields 

E{dr^,t) < Ce- Ao 'E(9 4 m ^,0). 

Summing on m yield the first estimate of the Corollary. The second follows from fl^Tl) . 
Sobolev imbedding, and the first. □ 
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3.2. Uniform A Priori Estimates. The existence of local-in-time solutions for the 
problem ( 1291) is proved as in [7j . In order to get global existence for the nonlinear problem 
( )29|) . it suffces to prove global a priori estimates for (129]) supplemented with initial data 



#(0,z) = hi{x), V t (Q,x) = h 2 {x). 
With C, Ao from Lemma HI choose T > so that 

a ■= Ce- XoT < 1. 
For t >T and integer k > 15, introduce 



(39) 



(40) 



[%*) 



\(Y,a)\ + \\(Y,a) 



where 



|(*,<r)| := sup 

re[0,t) 



E 



$^ e«||a;sr'*(^')lli»([„x|) + ei 



0<m<fc-6 \0<Km 



t/"V 



d£ r * 



and 



W,<t) 



sup 

0<T<t 



( 

E \\d l tdr^(r,.)\\ L2{[xo , L]) + E II#^*0v)||l>([.o,l]) 

. 0<Km, 0<l<k 

\0<m<k 

d k+l a 
+ sup ||d t fe+1 *(7v) - - M - T qi(;a)d x ^(r,-)\\ L 2 {[X0:L]) 
o<r<t a£^ x 



+ E (ii^r^(-^o)iu 2 [o, t ] + ii^9r^(-^)iiL 2[M )+ E 



0<l<m, 
0<m<k+l 



0<m<fc+l 



cPV 



d£ r 



Z»[0,t] 



with a positive constant A to be defined in (j52p below. 

In this subsection, when there are no specific indications to the contrary, we assume 
aij, bi, and g are functions of (x, \l/, V\l/, a, a), d\ and e\ are functions of (\l/, ^ t )- Define 



£($,£) := -f— (t,x ) + / ^+{(9 t $) 2 -a 11 (9 :E <l>) 2 }(t,x)da;. 

«lP+ </x 



(41) 



Furthermore, for any I G N and given \& and cr such that \\\(&, o~)\\\ < oo, define 

i ,i 

£,($,*) :=E)WM, and ' ^(*'*) := S-i(M+^(^-?iM**-j£,*) 



m=0 
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It is easy to see that if || |(\l/, a) | || < e for sufficiently small e > 0, then 
£($,*)(*) > C J {d t $) 2 + {d x $) 2 + ® 2 ){t,x)dx 
for some constant C > independent of £. 

Proposition 6. Assume that a' satisfies < T73]) . There exists an e > so £aa£ /or any 
< e < e 0; if (\1/, a) zs a smooth solution of the problem / Tlffij) and ( j3"9"|) satisfying 



then 



l°"o| + ll^i||tf fc + 2 + ll^lliffe+i < e 2 < Cq, and, 

111(^^)111 < c/2. 



f^crllll < e, 



(42) 



Proof: The proof has four steps. 

Step 1. Lower order energy estimate. Define 

@($,t) := / -u + (andi + a 01 )(dt®) 2 (T,x )dT+ / u + a 01 (d t ®) 2 (T, L)dr , 
Jo Jo 

and f^ m (<&,£) = X^/=n-^(^t^'^) - Taking the m-th (0 < m < k — 1) order derivative of 
the equation (1291) with respect to t, then multiplying the both sides by u + <9 i m+1 \I/ and 
integrating on Q =: [0, t] x [a;o,£], noticing that a o = 1 yields 



+ C\\W,a) 



e m 



m+2 



S m (^,r) + J2 



1=0 



d l a 



dt l 



dr + @ m (q,t) 



(43) 



for m — 0, 1, • • • , k — 1. 

Step 2. The highest order energy estimates. Take the k-th order derivative for the 
equation (12"9~1) with respect to £ to find, 

d fc+2 CT d'+V 



0<Kfc-l 



where 



& h {x, ¥, a) := X) °* j " E ^A'^ - E ^^V • 

l<Kfe I «j'=0 i=0 J 

In order to handle the term dtk+ 2 , study the equation for \1> = d£ \1/ — q±(x, a)^/ x -^-, 

Sf(x, *, a)^ = & k {x, *, o) + #(x, *, o), (44) 
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where 



#(*,*, a ) = J2 C] 



id 1 * 



a 



d k+l a 



0<l<k-l 



k-tf+2% l (Qi(x,a)^ x ) - ^^d t { qi (x,a)^ x ) 



d k a I \ { d k a 

- ~^d 2 t { qi {x,a)y x ) - I 2a m d t 8 x + a n d 2 x + s ^b i d i J l— qi ( x ,a)% 

Multiplying both sides of (I44p by u + d t ^ and integrating on Q yields 



+ C\\\{%a) 



where $} k is defined as follows 



At 

e 64 



fc+1 



£ h {V,T)+Y, 



1=0 



a 



dt l 



dT + @ k (V,t) 



d k a 



(45) 



@ k ($,t) = @ k ^($,t) + %(d?$ - qi (x,a)V x -- 7 ,t). 



dt k 



Adding the estimates (|4"3]) and (1451) yields 

&0M) + 40M) < 4(*,o) 

•t 



+ c-IIK^cr) 



At 

e 6 4 



fc+1 



4(*,r)+X) 



1=0 



d l o 



dt l 



dr + @ k {q,t) 



(46) 



Step 3. Boundedness of the energy. Differentiating the equation for the shock front 



a{t) = s/ 3 (*{t,x )) 



(47) 



with respect to t, yields 



fc+i 



1=0 



fc+i 



E W ^ C[\*{T)\* + y £\Bl*(T,X )f 



1=0 



Using ( l4"T|) again yields 

fe+i 



Thus 



z=o 



t k + l Al 

_at v^ da 

e 04 



d 1 



E^ ^ c[J2\d l Mr,xo)\ 2 + m\ 



.1=0 



n 



E 



Z=0 



dt l 



dr < C I / e"M^ (^,^^,^)^ + ^fc(^,t;^,a) 
Therefore, the energy estimate (j4"6"|) is equivalent to 

40M + 40M < 4(*,0) + C|||(*,<T)|||(# fc (*,r)+ f e-^£ k (^,T)dr). 
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If || | (*, a) 1 1| <e, then 



£ fc 0M) + ^OM < 2 4(^,0) < Ce 4 



This yields 



•1/2/ 



»V2/ 



(¥,<r)|| < C Bupfn^rJ + ^C*,*) < CV < 

,0<r« / I 



(4£ 



Step 4. Decay of the lower energy and the shock position. The basic idea to get the 
decay is to control the deviation of the solution \1/ of the nonlinear problem (j29p from that 
of the linearized problem ( 130]) (denoted by \I/). The contraction of the energy for \I/ will 
also yields the contraction of the energy for \1/. This gives the decay of the lower energy 
of \1/. The decay of the shock position is a consequence of the governing equation for the 
shock front and the decay of the lower energy. 

At time r = to, we can choose hi G H k and h 2 G H k ~ l so that there exists a solution 
* G C k - 1 - i ([t ,oo);H i ([x ,L])) of the linear problem (J3QJ satisfying *(t , •) = &i and 
^t^O) ') = ^2) an d ^ satisfies 

fc-i 2 
^£||^f(i ,-)l|zW] ^ C 111(^^)111 ( 49 ) 



Z=0 i=0 

for some uniform constant C, and 



£ fc _4(tt-tt,*o) < C|||(*,t7)|||f fc _ 4 (*,to). 
As that in Steps 1 and 2, energy estimates for the equation for $? — ^ gives 
4- 4 (* - *, *o + T) + 4- 4 (* - *, *o + T) - 4- 4 (* - *, t ) 



(50) 



<£ fc _4(tt-*,t )+C|||(tt,C7) 



to 



to+T fc ~3 

E 

Z=0 



dV 



dt J 



+ 4-4(* - * , r) dr 



;i/2 



;i/2 



- / 5 fc r 4 (*, r)C 4 (* - * , r)dr + (^_ 4 (*, *o + T) - <? fe - 4 (#, to)) 

'to 



+ (4-4(^-^,t + T) 



. 4 (*-*,t )) 



Using the contraction of energy for \l/ and noting that the deviation of \l/ and ^ at t 
is of higher order (cf. (1501 ). one has 



34 + 30a a /T ^ N 2 + 30ct a /T x 
£ fc _ 4 (^,£ + T) < -^£ fc _ 4 (#,£ ), 



64 



32 



(51) 
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if e is sufficiently small. As in Corollary |5l it follows from ( 151]) that 

£ k - 4 (%,t) + a 2 (t) < C(£ k ^,0) + a 2 (0))e- 2Xt , 
where with a from ( 140T) . 



Thus 



_ln( 1 + a )/2 ^ 2 + 30 a „ 

2T 17+15a 



fc_6 r jl 

a 1 a 



£ 



dtf 



< Ce 2 e- At . (53) 



+ ||*(*jVI|L-[.n,.L] 
(=0 

Combining ( 148]) and (153]) . one has (j42p . This finishes the proof of the Proposition [6j □ 



Once one has the Proposition [6] Theorem [2] follows from the standard continuation 
argument and local existence result. 
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